Experimental Investigation of the Evolution of Gaussian Quantum Discord in an 

Open System 
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Gaussian quantum discord is a measure of quantum correlations in Gaussian systems. Using 
Gaussian discord we quantify the quantum correlations of a bipartite entangled state and a separable 
two-mode mixture of coherent states. We experimentally analyze the effect of noise addition and 
dissipation on Gaussian discord and show that the former noise degrades the discord while the 
latter noise for some states leads to an increase of the discord. In particular, we experimentally 
demonstrate the near-death of discord by noisy evolution and its revival through dissipation. 
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Entanglement is undoubtedly a key resource in quan- 
tum information science and it has become synonymous 
with different nonclassical tasks such as quantum telepor- 
tation, dense coding and quantum computation How- 
ever, despite its obvious importance for such tasks, the 
exact need of entanglement in some other non-classical 
tasks has remained enigmatic. It has for example been 
shown that some quantum computational tasks based on 
a single qubit (the socalled DQC1 model) can be car- 
ried out by separable (that is, non-entangled) states that 
nonetheless carries non-classical correlations [2|-l4(. An- 
other quantum task requiring the use of quantum cor- 
relations but not entanglement is quantum key distribu- 
tion [5, 6]. This indicates that there exist genuine quan- 
tum correlations different from entanglement, which al- 
low for information processing that is intractable by clas- 
sical means. 

A measure of non-classicality that also captures the 
non-classical correlations of separable states is quantum 
discord as suggested in refs [7|, |8[ . The quest of under- 
standing the possible use of quantum discord for quan- 
tum tasks has fueled an explosion in the research of dis- 
cord spanning from the investigations of its evolution in 
noisy Markovian [9|-[l3| and non-Markovian [14|,ll5j chan- 
nels to its role in thermodynamics [l6[ and phase transi- 
tions in spin systems [17]. This measure has been given 
an operational interpretation in terms of the req uired re- 
source for enabling quantum state merging |l8l fl9j and 
local broadcasting |20j . 

Most work on quantum discord has been focused on 
finite dimensional systems including qubits and qudits 
but recently the discussion has been extended to infi- 
nite dimensional systems also known as continuous vari- 
able systems |2ll-[23| . This extension has mainly been re- 
stricted to the set of Gaussian states and Gaussian mea- 
surements [2l|, [22| , partly due to the associated reduced 
mathematical complexity of the problem and partly due 
to their immense importance for quantum information 
processing and their simplistic generation |24[ : Gaussian 
states can be efficiently prepared |25| and allow for secure 
quantum communication [6], near-optimal cloning |26j 
and teleportation [271 ] . 



In this Letter we experimentally characterize the Gaus- 
sian quantum discord of two-mode squeezed states and 
separable two-mode mixtures of coherent states. We in- 
vestigate its evolution in an open quantum system de- 
scribed by local Markovian decoherence in form of Gaus- 
sian noise addition corresponding to a classical noise 
channel and pure attenuation. 

The robustness of non-classical states against noise and 
loss is of high importance as real- world quantum informa- 
tion protocols will inevitably consist of noisy and lossy 
operations. Entanglement is known to be a fragile re- 
source which is difficult to generate and the optimal use 
of it would require complex purification and distillation 
protocols. On the contrary, it has been shown theo- 
retically that quantum discord of some states is robust 
against Markovian decoherence, and in fact discord can 
increase through dissipation both for discrete [9l4l l\ and 
for continuous [12[ variable systems. By interrogating 
two-mode mixtures of coherent states in an open system 
we show experimentally that Gaussian quantum discord 
is reduced in a classically noisy channel and, counterin- 
tuitively, increased in a simple dissipative channel. 

Figure [1] shows the experimental setup that was used 
to generate quantum correlated Gaussian states. A pair 
of optical parametric amplifiers (OPAs) based on type 
I quasi-phase-matched periodically poled KTP cystals 
placed inside bow-tie shaped cavities were used to gener- 
ate two independent amplitude squeezed beams at 1064 
nm. The OPAs were seeded with dim laser beams at 
1064 nm to facilitate the locking of the cavities and sev- 
eral phases of the experiment. Both OPA outputs had 3.2 
±0.2 dB of squeezing and 6.7 ±0.2 dB of anti-squeezing 
which was measured with a homodyne detector with a 
total efficiency of 85 ± 5% and electronic noise contri- 
bution of —20 dB relative to shotnoise. The measure- 
ments were performed at the sideband frequency of 4.9 
MHz with a bandwidth of 90 kHz. In order to generate 
correlations in addition to those produced in the OPAs, 
one of the seed beams was symmetrically modulated with 
two electro-optic modulators - an amplitude and a phase 
modulator - that were driven by independent electronic 
noise generators. We interfered the two squeezed beams 
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on a symmetric beam splitter with a relative and locked 
phase of tt/2 to produce a pair of quadrature entangled 
beams. One of the output beams propagated through a 
dissipative channel which was implemented by a beam 
splitter with a variable transmittivity. This setup can 
hence generate and characterize a variety of two-mode 
Gaussian states. 



State Preparation 



OPAl 



SH Pump 




FIG. 1: Experimental setup. See main text for description. 
Abbreviations: IR: infrared (1064 nm), SH: second harmonic 
(532 nm), PM and AM: phase and an amplitude electro opti- 
cal modulator, OPA: optical parametric amplifier, BS: beam- 
splitter, LO: local oscillator. DAQ/SP data acquisition/ sig- 
nal processing 



Due to the Gaussian nature of the generated two-mode 
states, they are fully characterised by their covariance 
matrices which are easily estimated from the homodyne 
measurements of the amplitude quadrature X, and phase 
quadrature P. The covariance matrix for the state p AB 
written in the standard form is 
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where 7 = Diag [cov(x A , xb), cov(pa,Pb)], a = 
Diag [var(x^), var(p^)] and f3 — Diag [var(x j e), var(ps)]. 

Gaussian quantum correlations beyond entanglement 
are captured by the measure of Gaussian discord which 
we formalise in the following [2l|, [22|. In a bipartite 
system, the total amount of correlations (classical and 
quantum) is given by the von Neumann mutual infor- 
mation I(pab) = S(Pa) + S(Pb) - S(p AB ) where S(p) is 
the von Neumann entropy and p A(B) is the reduced den- 
sity matrix of the A (B) subsystem. Another measure of 
mutual information that only quantifies the amount of 
classical correlations extractable by a Gaussian measure- 
ment is J A (p AB ) = S(p A ) - inf aM S(p AWM ) where a M is 
the covariance matrix of the measurement on mode B. 
As it only captures the classical correlations, the differ- 
ence, V A = I(pab) — Ja(Pab), is a measure of Gaussian 
quantum correlation that is coined Gaussian quantum 
discord. An explicit expression for this discord has been 
found IH: 
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FIG. 2: a) Mutual information for the two- mode squeezed 
state with the circles representing the / information (red) and 
the J information (blue). The curves correpond to theoreti- 
cal curves fitted to the first data point, b) Quantum discord 
and logaritmic negativity (inset) for the two-mode state as 
a function of attenuation. The experimental data (circles) 
show the discord and logaritmic negativity calculated from 
the measured covariance matrices and the solid lines are the- 
oretical fits to the first measurement point (corresponding to 
a channel attenuation of 0%). 



V(a AB ) = f( 02) - f(iz-) - f(y+) + f(VE^) 
where f(x) = (^l)log [2±±] - (^)log [^] and 
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where a) applies if (J 4 - hh) 2 < ll(h + l)(ii + h) 
and b) applies otherwise. I\ = det[a], I2 = det\j3], Is = 
det[y], I4 = det[cr AB ] are the symplectic invariants and 
v\ = \{5 ± VS 2 - 4/4), with S = h + h + 2/3, are the 
symplectic eigenvalues. 
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Using the experimentally obtained covariance matrix 
of the two-mode squeezed state, we calculate the Gaus- 
sian discord for different attenuations of mode B. The 
results for the mutual informations are shown in Fig- 
ure^ and the discord in Figure^. As usually expected 
for quantum correlations in a dissipative channel, we see 
a monotonic decrease of the Gaussian discord with in- 
creasing dissipation. A similar behavior is expected and 
observed for the entanglement of the state as shown in 
the inset of Figure [2] where the logarithmic negativity is 
plotted against attenuation [28|. We thus conclude that 
for this particular two-mode squeezed state, dissipation 
has a degrading influence on both entanglement and dis- 
cord. However, it is well-known that entanglement is not 
the only resource containing Gaussian discord. A two- 
mode mixed separable state may also contain Gaussian 
discord, and as we shall see in the following, this discord 
is much more resilient to dissipation than the entangle- 
ment based discord. 

In the experiment we generate a separable mixed state 
by enabling the electronic noise generators and disabling 
the OPA's (see Figure [TJ. This corresponds to the split- 
ting of thermal state on a symmetric beam splitter and 
thus the generation of a two-mode mixture of coherent 
states with correlations between the amplitude quadra- 
tures and the phase quadratures of the two modes. For 
this state, we measure the covariance matrix for different 
modulation depths of the modulators and subsequently 
calculate the discord with the results shown in Figure [3] 
(green circles). The theoretically expected behavior of 
the discord for the mixed state is represented by the 
solid red line which is monotonically increasing, even- 
tually saturating for large modulations. We see that the 
experimentally obtained values of the discord decrease 
for very large modulations thus deviating from the the- 
oretically predicted behavior. This discrepancy is sim- 
ply due to finite balancing between the two homodyne 
detectors which was measured to have a common mode 
rejection ratio of 27 dB. Such an imperfection of our de- 
tection system leads to the measurement of uncorrelated 
events which directly simulates a classically noisy chan- 
nel where the amount of noise scales with the modulation 
depth. Since we are interested in simulating the effect of 
classical noise on the discord, we decided to reduce the 
common mode rejection ratio of the two detectors to 15 
dB. The results for the discord are shown by the red 
circles in Figure |3Jl, and the constituents - the mutual 
informations / and J - are presented in Figure [3]3. In 
a classically noisy channel we thus observe the transi- 
tion from a quantum regime with non-zero discord to a 
classical regime with near zero discord. 

In the next step we investigate the effect of dissipa- 
tion on mode B. Experimentally, dissipation is simulated 
by a simple beam splitter placed in mode B as shown in 
Figure [TJ The Gaussian discord as a function of channel 
attenuation is estimated through measurements of the 
covariance matrix, and the results are shown in Figure |4] 
for different amounts of modulation (partially resulting 



a) 07 



0.6 



0.5 



0.3 



0.2 



O Experiment 15 dB CMR 
O Experiment 27 dB CMR 

Ideal theory 

- -Theory for 15 dB CMR 
Theory for 27dB CMR/ 



7 n 



0, 



b) 



5 10 
Modulation depth [db] 



15 



20 



O Experiment 1 5 dB CMR / 

° Experiment 1 5 dB CMR J 

Ideal theory / 

Ideal theory J 

_ 15 dB CMR theory/ 

co 4 - - - 1 5 dB CMR theory J 




5 10 15 

Modulation depth [db] 

FIG. 3: a) Quantum discord as a function of the modula- 
tion depth. The circles represent the experimentally obtained 
data for discord with 27 dB (green) and 15 dB (blue) common- 
mode-rejection (CMR) of the homodyne detectors, b) Mutual 
informations / and J as a function of the modulation depth. 
The circles represent the experimentally obtained data. Solid 
lines result from theory with perfect detectors (corresponding 
to a noise-free channel) whereas the dashed and dot-dashed 
curves are the results of the theory including CMR and elec- 
tronic noise. The error bars are given by the statistical un- 
certainties. 



in classically noisy channels). From these measurements 
we see, interestingly, that the discord increases as a func- 
tion of local and Markovian dissipation. We do not only 
observe an extreme robustness of discord associated with 
the separable state but also a clear rise of quantum cor- 
relations during dissipative dynamics. This is in stark 
contrast to entanglement which cannot increase under 
any local (unitary or non-unitary) transformations (29j . 
Ultimately, it is possible to observe the complete death 
of discord through classical noise addition and its sub- 
sequent revival through dissipation as partially demon- 
strated in Figure [3] (near death) and Figure 2] (revival) . 
We interestingly note that the discord of the classically 
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noisy state investigated above can be revived after death 
even if the state at A has been measured, stored in a 
classical memory and subsequently recreated. This is 
also clear from the fact that a classically noisy channel is 
theoretically identical to a measurement map [20]. 
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FIG. 4: The discord of four different mixed states from Fig- 
ure [3] with different modulation depths as a function of atten- 
uation of mode B. The error-bars are given by the statistical 
uncertainty. The solid lines are theoretical estimations, fitted 
to the first data point associated with zero attenuation. 



The revival of the discord as a result of dissipation can 
be understood from two effects. The first one is the at- 
tenuation of the uncorrelated classical noise in mode B, 
which is responsible for the near death of the discord in 
Figure [3] This effect is hence the main contributor to 
the revival of the discord when uncorrelated noise above 
shotnoise is present. The second effect that leads to a 
further increase in the discord results from the relatively 
higher amplitude of mode A. This makes the inevitable 
noise contribution of a measurement on mode B more sig- 
nificant and thereby increase the difference between the 
information obtainable by quantum and classical means. 

Although the Gaussian discord of our experimentally 
obtained separable states is more resilient to losses than 
that of the entangled states, it is important to point out 
that for the same amount of energy of the two types of 
states, the entangled state will for any amount of Marko- 
vian loss exhibit the highest discord. This is illustrated 
in Figure [5] where the discord of pure entangled and sep- 
arable states with mean photon numbers between 1 and 
100 (1 and 10 photons corresponding to 5.7 dB and 13.4 



dB of two-mode squeezing) is plotted against attenua- 
tion (in dB). Needless to say, such pure squeezed states 
are much more difficult to produce than the more robust 
two- mode mixed coherent state with 100 photons (green 
curve in Figure [5j) which is straightforwardly produced 
with electro-optical modulators. 

Finally, we note that the discord of a combination of 
entangled states and separable states might be the opti- 
mal resource for some applications. E.g. we have recently 
shown that by exploiting the correlations of combined 
two-mode squeezed states and coherent mixed states the 
performance of a continuous variable quantum key dis- 
tribution system can be improved in terms of increased 
robustness to channel noise and dissipation (30| . 
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FIG. 5: Theoretical simulations of the evolution of Gaussian 
discord in an attenuating channel for pure two-mode squeezed 
states and two-mode mixtures of coherent states. 



In conclusion, we have experimentally generated two 
different states containing Gaussian quantum discord: 
A two-mode entangled state and a two-mode separa- 
ble state. The discord associated with the entangled 
state was shown to degrade under dissipative dynamics 
whereas the discord of the separable state was exhibit- 
ing the opposite behavior; it was an increasing function 
of dissipation until a certain point after which it rapidly 
vanished. Furthermore, we show that discord of separa- 
ble states is reduced if uncorrelated noise is added to the 
state but that it can be revived through dissipation. 

We acknowledge the financial support from the Danish 
Council for Technology and Production Sciences (no. 10- 
093584 and no. 10-081599). 
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